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Abstract 

In this paper we theoretically show that the Raman scattering by a core-shell micron or submi- 
cron particle with epsilon-near-zero metamaterial shell and silica spherical or cylindrical core can 
combine useful features of cavity-enhanced and surface-enhanced Raman scattering. The cavity 
resonance together with the plasmon resonance lead to the giant enhancement of the field inside 
the metashell which is performed as a layer of silver or gold nanoparticles and is penetrable for 
molecules to be detected. This approach results in the significant increase of both effective vol- 
ume in which molecules are affected by enhanced electric field and Raman gain averaged over this 
volume. 

PACS numbers: 78.45. Ih, 78.30.Er, 78.30.Fs, 78.67.Bf, 78.55.Et, 73.20.Mf, 85.352p, 05.45.-a, 07.79.Fc, 
42.55.Sa 
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The role of the surface-enhanced Raman scattering (SERS) \^ in the modern sensing, 
especially in molecular detection, is huge. The mechanism of SERS is related with plasmonic 
nanoparticles resulting (in conventional SERS schemes) from roughening the silver interface. 
The nanoparticles offer the resonant enhancement of the local field acting on molecules 
(located near them in a liquid or gas host medium) compared to the incident wave field. This 
enhancement corresponds to the proportional increase of the molecule dipole moments at the 
excitation frequency 0e and at Raman frequencies 0r (e.g. Q). This effect is complemented 
by the similar increase of the radiation of the pair molecule plus nearest nanoparticle at 
the Raman frequency (e.g. . The plasmon resonance experienced by a plasmonic 

nanoparticle (sphere or ellipsoid) is rather wide-band. Therefore, usually, one of Raman 
frequencies radiated by molecules or both of them together with the excitation frequency 
lie within the nanoparticles resonance band. Practically, what is detected in SERS is not 
the radiation of molecules at frequency 0r but the re-radiation of resonant nanoparticles 
excited by molecules at 0r. For the amplitude of the field scattered by the pair molecule 
plus nanoparticle enhanced due to the presence of the nanoparticle the coefficient K{0e) 
expressing the local field enhancement is multiplied by the coefficient x{0r) which equals 
to the ratio of dipole moments of a resonant nanoparticle to that of a molecule. In the 
classical theory 2, s], 4| it is shown that k{0) = x(0)- When the Raman shift is small i.e. 
0_R ~ 0e ~ 0av = 0e/2 + 0i?/2 this rcsult corrcspouds to the SERS amplitude electromagnetic 
gain Ginax = ft:max(0e)Xmax)(0i?) ~ '«max(0av) of the ordcr IQ^^'^-^.s) ^^le SERS intensity 

electromagnetic GInvAx = ^max(0 

av) ~ 10^"^ ^-^ [2; [sj . Here the subscript corresponds to 
the optimal location of a molecule with respect to the nanoparticle and the incident wave 
vector because only the radial (with respect to the nanoparticle center) field component 
is enhanced, for polar and azimuthal components of local field we have Gmax < 1- The 
interval Kjnax(0av) = 3 . . . 10 corresponds to the characteristic interval of distances between 
the molecule and the nanoparticle surface (usually 2 — 3 nm) 0, Q]. In random arrays of 
nanoparticles their electromagnetic interaction changes the electromagnetic SERS gain Gmax 
not significantly. 

Some other mechanisms (chemical adsorption of the metal interface that modifies the 
molecule polarizability [5] and some quantum effects [6]) bring additional increase of the 
molecule dipole moment at Raman frequencies. As a result, SERS for a molecule located 
at same distances from the surface of a spherical nanoparticle corresponds to the total 
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amplitude gain G 



tot 



;i - 3) ■ 10^ or total intensity gain GP^l^ ~ lO^*'-^) 



Qfl. 



Here the 



superscript means that in coefficients Gmax('^av) and G/^*j^(0av) not only electromagnetic 
mechanisms of SERS are taken into account. 

Specially grown nanowires or nanodisks instead of spheres or ellipsoids particles improve 
the result for the maximal electromagnetic gain 10^ . . . 10^ times for points near the wires 
ends or disk edges (e.g. [71]). Tremendous electromagnetic gain G/max ~ 10^^ can be obtained 
for molecules located within small (0.5 — 2 nm sized) gaps between paired nanoparticles 
(nanospheres, nanorods, bow-tie arms) as well as in similar gaps between almost touching 

u\~\ 

particles of plasmonic nanoclusters [9|, 110[ . That allowed one to apply SERS for detecting 
separate molecules jol. 

A drawback of these exciting variants of SERS is the extreme locality of the huge field 



enhancement [111 . Il2l |. Thus, this method practically requires to direct a molecule under 
detecting to a selected point. In [8| the practical importance of the highest possible averaged 
electromagnetic SERS gain G/av was stressed. The averaging should be done over the 
effective domain occupied by metal nanoparticles or nanocorrugations (minus the volume 
occupied by the metal) j^. In values GI.^^ ~ 10'' were theoretically engineered in a regular 
array of exactly touching cylindrical nanocorrugations on the silver half-space. Since [&] the 
progress in the design of high G/av has been significant. The mechanism of the high G/av 
is related to the electromagnetic interaction in regular arra^ of plasmonic nanoelements. 



In regular arrays values of G/av can attain G/a 

HE? 



2 ■ 10^^ 



13| that is accompanied by the 



maximal gain Gl-a 
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14-15 



llq at the crevices between the corrugations . However, 
these amazing values of the averaged gain require the extreme precision of nanofabrication. 
Statistical deviations in the array geometry of the order of one Angstrom lead to its dramatic 
decrease from G/av ~ 10" to G/av ~ lO^^"^) [isl. 

The idea of the present paper is to show the way to very high values of G/av using random 
arrays of metal nanoparticles on the dielectric core to be fabricated by the self-assembly. 
Very high G/av theoretically results from the combination of SERS with cavity-enhanced 
Raman scattering. The last one is an important direction of the modern literature. The 



strong field enhancement holds around the points corresponding to the w 
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23| that 



(WG) modes maxima inside the optical microcavity 
also results in the high Raman gain (defined for these structures in a different way than in 
SERS). Definitely, this method is not applicable for detecting the separate molecules since 
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the field of WG modes is concentrated inside tlie solid cavity. The cavity-enhanced Raman 
scattering effect is practically thought as promising for microlasers and is considered usually 
as the stimulated Raman scattering (SRS) effect |24|. In SRS the incident light is used for 
pumping the WG states of the microcavity which are excited at Raman frequencies of the 
cavity material. 

Arrays of touching microcavities were suggested and studied in work 25(] with the purpose 
to enhance the local field around the contact points. This approach combines some features 
of SERS (sensing of molecules in the host medium) and cavity-enhanced Raman scattering 
(WG modes at 0r). However the result for G/av (defined in this case the same way as in 



SERS) is as modest as GI^v = 7 25|. Notice that the electromagnetic coupling of an optical 
microcavity (their typical radiation quality is as high as Q ~ lO'''^"^"-') to outer wave fields is 

21| 



negligibly weak. Thus, the simple pumping by incident waves is not efficient [18|, [19 
and in cavity-enhanced Raman scattering schemes one uses special coupling elements 
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prisms 



25|. 



or waveguides) with wave leakage between them and cavities 

The goal of the present paper is to show that using a layer of silver nanoparticles (e.g. 
spheres) randomly distributed on the silica core of submicron or micron size as shown in 
Fig. [1] (a) one can significantly improve G/av compared to the conventional SERS based on 
the same nanoarray on the planar substrate. Moreover, using such core-shell clusters one 
increases the effective volume where the field is enhanced. Such "templated nanoparticles" 
(in the terminology of works 29, 30|) can be prepared with existing technologies 29. 30. Isi] 



which allow one to obtain also the cores with high precision of the diameter [28||. The 
technology of preparing dense one- layer and even multilayer arrays of these "templated 
nanoparticles" (or "templated microparticles" ) on a dielectric substrate is described in [29|, 



30|. Compared to a a conventional SERS scheme with plasmonic nanoarray covering a planar 



substrate the use of spherical "templated particles" gives the gain in the total effective 
volume of the nanoarray approximately equal to 27r. The use o f cy lindrical cores gives 
the volume gain close to tt. This fact was noticed in works 29|, l30|] where the SERS in 
such core-shell nanoclusters was experimentally studied. The experimentally demonstrated 
averaged (over the metashell volume) Raman gain was nearly the same as in conventional 
SERS with the same nanoarray on the planar substrate. The best known result G/*"* ~ 10^ 



has been obtained in 



3l|. Below we predict much better results for optimized parameters 



of " templated particles" . 
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In the present paper we consider "templated particles" with spherical or cylindrical ge- 
ometry. The metashell couples the cavity with the host medium and the WG modes can be 
efficiently excited by an incident plane wave. Our design goal is to engineer the cavity reso- 
nances so that the WG modes would be concentrated inside the metashell. The main design 
parameter is the core radius a. Though the metashell is the same plasmonic nanoarrays 
which is used in conventional SERS (on planar substrates), the averaged gain significantly 
improves due to the WG resonance. 

The metashell can be presented with high accuracy as a layer of an effective metamaterial 



shown in Fig. [T] (b) [26|, l27|, |3^ . In works [26|] the homogenization model of the spherical 
metamaterial samples with radius a = 20 — 50 nrn (diameters of nanoparticles c? = 3 — 6 
nm) was validated by additional calculations. In 27| the same was done for a metashell 
on a silica core. In 3^ the high accuracy of the homogenization model for the metashell 
was confirmed by measured optical spectra where diameters of nanoparticles was either 
c? = 15 or c? = 20 nm and the core radius was a = 135 nm. These results allow us to use 



analytical calculations basec 
a concentric layered sphere 



on the known solutions of two plane-wave diffraction problems: 



3^1 



for the fie^ 
found in 



and a concentric layered cylinder 3J|. Explicit expressions 
d in the inner (core), intermediate (shell) and outer (host medium) regions can be 
35| for a sphere and in 36| for a cylinder (these papers were used for testing the 



Matlab codes). Both Maxwell Garnett (as in |27|]) and Bruggeman (as in 32]) models were 
used for the metashell of silver nanospheres. The chosen design parameters of the metashell 
(the nanosphere diameter d = 6 nm and the averaged gap between nanospheres 6 = 1 nm) 
correspond to the intervals where the Maxwell Garnett and Bruggeman models give the 
same values of the metashell effective permittivity e. The host medium permittivity was 



37( 1 and that of silver 



assumed to be free space. The permittivity of silicon was taken from 
from [38|. 

Calculations of the electric field showed that the WG modes excited by the incident plane 
wave concentrate inside the metashell at the blue edge of the collective plasmon resonance 
of the metashell (where its permittivity is close to 0). 

First, let us prove that the formula Gav = '^av(0e)Xav)(0R) 



'«max(0av) Or Gl 



av V i^av 7 



still holds for "templated particles" when 0e ^ 0r 



,. Let a molecule be located 



at a point A inside the metashell of the spherical core. The radial component of its dipole 
moment at the Raman frequency is denoted in Fig. [T] (b) as pi. Two other components of the 
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FIG. 1: (Color online) (a) - A silica micron- or submicron-sized sphere (or a cylinder) covered 
with metal nanoparticles. (b) - The same structure presented as a core-shell cavity. A point dipole 
pi located in the metashell (the radial component of the molecule dipole moment) creates at the 
observation point the which equals to the radial component of the field produced at the center of 
the dipole pi by an auxiliary dipole p2 = Pi placed at the observation point. Eq is the amplitude 
of the field of the same auxiliary dipole in the absence of the core-shell cavity. 

dipole moment are not significant for SERS in spherical "templated particles". This dipole 
creates at an arbitrary chosen observation point B located in the far zone the 6'-polarized 
field of complex amplitude E{B) that we can present in the form E{B) = x{A)pF{R). Here 
pF{R) {F is the known function) expresses the field produced at point B in absence of the 
"templated particle" by the y-oriented dipole p = pi located at the same distance R from 
B as the "templated particle" center (see Fig. [1]). From reciprocity E{B) equals to the 
radial component of the field E'{A) produced at point A by an auxiliary y-oriented dipole 
P2 = Pi = P located at B. Since the dipole p2 is located very far from the "templated particle" 
the incident field from which the radial field E{A) = El{A) results can be approximated as 
the field of a plane wave with amplitude Eq = F{R)p at the sphere center. By definition 
the local field amplitude enhancement is equal to k,{A) = \E{A)/Eq\, where as we have seen 
E{A) results from the plane wave refraction. Therefore the amplitude radiation enhancement 
due to the presence of the "templated particle" defined as = \E{B)/pF{R)\ (i.e. with 

respect to the same molecule located at the same distance R from the observation point as the 
"templated particle" center) is equal to |x(A)| = \E{B)/pF{R)\ = \E{A)/Eo\ = k{A). Since 
this result holds for an arbitrary point A, the averaged radiation and local field enhancements 
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FIG. 2: (Color online) (a) - Wide-band frequency plot of the averaged local field intensity gain G 
(in dB) for spherical core with radius a = 323 nm and of the complex permittivity of the metashell. 
|(b)| - Same as (a) in a narrow band for three values of the core radius a = 322 nm (solid line), 
a = 323 nm (dashed line), a = 324 nm (dotted line). 



are also equivalent. So, the enhancement of the molecule radiation at the Raman frequency 
in presence of the spherical "templated particle" is equal to the local field enhancement (with 
respect to an incident plane wave) in presence of the same nanoparticle in the same place. 
The same speculation can be done for other field components, for a cylindrical "templated 
particle" , and for any other finite structure modifying the local field. 

To find Gav we have to average the intensity |-E'(y4)p over the metashell volume (for 
spherical and cylindrical cores it is a simple numerical integration of analytical expressions). 
Here E{A) = K(0av, A) is the complex amplitude of the radial component of electric field 
when the "templated particle" is impinged by a plane wave of unit amplitude. The intensity 
Raman gain is equal to G/av = Gl^. 

The result for Gav (in dB) versus light wavelength for the spherical "templated particle" 
with the core of radius a = 323 nm is shown in Fig. [2] (a) together with the complex 
permittivity of the metashell. At the collective plasmon resonance of the metashell, the 
losses are very high and G^v ^ 1. The "useful" WG resonance (when the field concentrates 
inside the shell) corresponds to the azimuthal and polar numbers of the WG mode L = 1 
and = 9, respectively. It holds at A = 366 nm where e = —0.1 +i0.16 and Gav = 5784. It 
corresponds to the averaged intensity Raman gain G/av ~ 4 ■ 10^ and to the cavity optical 
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quality Q ^ 5 ■ 10^. The presence of the metashell broadens the band of the WG resonator. 
Cores with a = 322 and a = 324 nm still correspond to Glav ~ 3 ■ 10^, as we can see from 
Fig. [2] (b), where (logj^gGav is shown in more details for three values of the Si sphere radius. 
This result shows the possible tolerances in the cavity fabrication. 

It had been expected that the higher-order WG resonances (when a > 1 fim) should 
have strongly increased the gain since the optical quality of cavities usually grows along 
with the resonance order. However, for "templated spheres" it is not so. The next "useful" 
WG resonance corresponds to a = 611 nm (L = 1, = 18), however G/av ~ 3 does not 
improve since the WG resonance shifts to lower frequencies where lm{e) increases and the 
losses reduce the gain. 

For cylinders this blue shift is absent and the gain grows versus the order of the " useful" 
resonance. The case shown in Fig. [3] (a), i.e. the normal incidence of an axially polarized 
plane wave to an infinitely long cylinder of radius a covered with the same metashell was 
studied. There is only one component of the electric field E = E^. In Fig. [3] (b) the 
distribution of the field intensity over the structure with a = 337.5 nm at A = 353.5 nm (the 
mode azimuthal number L = 9) is shown. The localization of the field inside the metashell 
is clearly seen. One can see that at this frequency Gmax = I -^^ I max exceeds 500. However, 
the numerical averaging gives Gav = 202 i.e. G/av ~ 4 ■ 10^. 

Fig. m (a) shows the frequency dependence of Gav for three values of the Si cylinder 
radius a = 615.5, 616, 616.5 nm corresponding to the WG modes also concentrated within 
the shell (the resonance order L = 18). Then the result is much better Gav = 2500 . . .2700 
(G/av ~ 10'') within the band A = Aav ± 1 nm, where Aav = 352 nm. In Fig. H] (b) the results 
are shown for an abstract metashell with dispersionless permittivity ^const = —0.08 + ^0.156 
which is equal to £:(Aav)- It is clear that the dispersion of the metashell permittivity has no 
impact but for the highest gain one has to adjust e (in the region of \e\ <^ 1). The same 
structure at same frequency gives Gav = 3 • 10^ for £const = —0.01 + iO.Ol. 

Further increase of the cylinder radius allows us to obtain the higher gain. The result 
corresponding to the WG resonance of the order L = 32 is Gav(A = 353nm) ^ 10^. In this 
case the band of the WG resonance becomes as narrow as 10"'^ nm and the optical quality 
is of the order Q ~ 10''. Notice, that very high values of G/av correspond to very narrow 
resonance band. Then the excitation frequency and the Raman frequency can belong to 
this resonance band only if the Raman shift is very small. Otherwise the approximation 
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FIG. 3: (Color online) (a) - Cylindrical silica core with the same metashell impinged by incident 



plane wave, (b) - Electric intensity distribution in the structure with a = 337.5 nm at A = 353 nm 
and Gav = 202. 



0^ ^ 0^ is not anymore valid, and the result |x(0k)| = /t(0e) does not hold. The 

"template particles" with a > 1 /im have at the "useful" WG resonance the optical quality 
of the order Q > 10^ and huge Raman gain G'/av(Aav) ~ 10^^. However this gain can be 
applied only for molecules with very small Raman shift (less 1 cm~^ for wavenumbers). 
Therefore we do not consider so big "templated particles" in more details. 

To conclude: it is theoretically demonstrated that the WG resonances in core-shell "tem- 
plated particles" whose effective shell permittivity is close to zero are very promising for 
SERS. These structures combine useful features of SERS and cavity-enhanced Raman scat- 
tering. Exact analytical calculations were also done (for both spheres and cylinders) covered 
with a solid silver shell of same thickness d = Q nm. The WG resonances at which the field 
of WG modes is strongly localized (concentrated at the inner interface between the core and 
the shell) were found. These resonances demonstrated tremendous Gav that attains even for 
submicron a values Gav ~ 10^ (G/av ~ 10^^). This effect can be applied in prospective SRS 
schemes. 
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